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Abstract 

The perturbative contribution in the average plaquette is subtracted using Borel summation and the remnant of the 
plaquette is shown to scale as a dim-4 condensate. A critical review is presented of the renormalon subtraction scheme 
that claimed a dim-2 condensate. The extracted gluon condensate is compared with the latest result employing high 
order (35-loop) calculation in the stochastic perturbation theory. 
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1. Introduction 

Extracting the gluon condensate from the average pla¬ 
quette is an old problem of lattice gauge theory. The 
operator product expansion (OPE) for the plaquette in 
pure SU(3) Yang-Mills theory is given by 

PiP) = (l-iTrUn) 

= y ^ + ^Z06)<-GG)a4 + 6>(fl6), (1) 
36 n 

n=l ‘ 

where p is the lattice coupling and a the lattice spacing. 

To extract the gluon condensate the perturbative con¬ 
tribution, which dominates the plaquette, must be sub¬ 
tracted accurately, to better than one part in lO'*. The 
perturbative series is expected to be an asymptotic se¬ 
ries and at high orders the coefficients to be dominated 
by renormalon-caused large order behavior. However, 
the perturbative coefficients computed using stochastic 
perturbation theory 111] to 10-loop order does not display 
a renormalon behavior but a power law that grows much 
faster than the expected renormalon-caused large order 
behavior. 

To subtract the perturbative contribution Burgio et al. 
introduced a continuum scheme in which the renor¬ 
malon contribution to the perturbative part of the pla¬ 
quette is calculated by matching the computed high or¬ 
der coefficients with the renormalon-caused large order 
behavior mapped to the lattice scheme. This scheme of 
subtracting perturbative contribution resulted in a power 
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correction that scales, surprisingly, as a dim-2 conden¬ 
sate. This is in contradiction with the OPE ([Til that de¬ 
mands the leading power correction be of a dim-4 con¬ 
densate. 

In this presentation based on |[3l we give a criti¬ 
cal review of the above renormalon subtraction scheme 
and show that the continuum scheme in which the 
renormalon was subtracted is far from a renormalon- 
dominated scheme and it fails as well the self- 
consistency test. 

We then introduce a renormalon subtraction scheme 
based on the bilocal expansion of the Borel transform in 
continuum scheme, and show that the plaquette data less 
the Borel-summed perturbative contribution scales as a 
dim-4 condensate. The extracted gluon condensate is 
then compared with the result from the recent high order 
(35-loop) calculation of the plaquette using stochastic 
perturbation theory. 

2. Critical review of existing renormalon subtrac¬ 
tion scheme 

The renormalon subtraction scheme of Burgio et at. ||2l 
writes essentially the average plaquette as 

P(P) = P‘™(/3,) + dPiPc) + PnpCjS) , (2) 

where 

SJ 

with Pc denoting the coupling in the continuum scheme 
defined by 

Pc^p-n-'j (4) 
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Here the plaquette is divided into the renormalon contri¬ 
bution F™ and the rest of the perturbative contribution 
i5F, and nonperturbative correction Pnp- The asymp¬ 
totically divergent behavior of the perturbative series is 
contained in F™, and dF denotes the rest that can be 
expressed as a convergent series. 

Now define with 


N 

F^^^(/3) = F(J3) - - CDp-'’, (6) 

n=l 


where C™ denotes the perturbative coefficients in the 
lattice scheme of Note that, by definition, P[^^ is 
free of perturbative coefficients to order N. The con¬ 
stants ri, r 2 that define the continuum scheme and the 
normalization constant M are to be determined so that 
CJ,®" converges to c„ at large orders. In the continuum 
scheme given by 


coefficients between the lattice scheme and the contin¬ 
uum scheme of O at the orders in consideration are still 
very sensitive on the low order coefficients, which were 
ignored in the matching procedure of |12l . It is thus obvi¬ 
ous that the continuum scheme of O cannot be the right 
scheme where renormalon can be subtracted reliably. 


Table 1: The peiturbative coefficients of the average plaquette in the continuum 
scheme. 
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Checking the internal consistency of the subtraction 
scheme also shows an underlying problem. The nonper¬ 
turbative term in (|2|l can be written using (|6]l as 


PNpOS)=Pjfp> 


08)- 


6F(fi,)-Yicn-CTW 


n=l 


( 8 ) 


For P[^' to approximate the power correction 


ri=3.1, r2=2.0, (7) 

and with an appropriate value for N it was observed that 
CJ,®" converge to c„ at the orders computed in stochas¬ 
tic perturbation theory. Because the last term in (|6]) is 
a convergent series P[^^ is well-defined at » cxi, and 
this is precisely the quantity that was assumed to rep- 
resent the power correction, and it was P^^p that was 
shown to scale as a dim-2 condensate. 

In this procedure the large order behaviors in the lat¬ 
tice scheme and the continuum scheme are matched 
using only the high order coefficients. However, this 
matching of large order behaviors would only work 
when the perturbative coefficients in both schemes are 
already in asymptotic regime and follow the pattern of 
renormalon-caused large order behavior. But, because 
the computed coefficients in the lattice scheme are far 
from being in the asymptotic regime and do not fol¬ 
low the renormalon pattern the matching cannot be per¬ 
formed. Therefore, the conclusion of a dimension-2 
condensate based on this matching should be reexam¬ 
ined. 

One way to check the consistency of the above pro¬ 
cedure is to map the computed coefficients in the lattice 
scheme to the continuum scheme and see if the mapped 
coefficients follow the pattern of a renormalon behavior. 
As can be seen in Table [1] however, the coefficients in 
the continuum scheme are sign-alternating, instead of a 
factorially growing pattern expected from a renormalon 
behavior. It shows that the mapping of the perturbative 


N 

6F(/3c) - ^(c„ - C'®“)^-" 

n=\ 


« <pC8) 


(9) 


must be satisfied. Because 6F(J3c) is by definition a con¬ 
vergent quantity it can be written in a series expansion 


SF(p,) = YDnl3c''^ ( 10 ) 

n=l 


where D„ can be computed to the order c„ are known, 
and (|9]l can be written approximately as 

ii:tig«y8-"-i:ii(c„-cr);8-"i 

<pV) ■ ^ ^ 

In the continuum scheme of (jT), and at A = 8 and 
P - 6.0,6.2 and 6.4, for example, the ratios are 69,59 
and 42, respectively, which is a severe violation of the 
consistency condition. This again confirms that the con¬ 
tinuum scheme of Q cannot be a proper scheme for 
renormalon subtraction. 


3. Renormalon subtraction by Borel summation 

Clearly, the perturbative contribution in the plaque¬ 
tte cannot be subtracted by matching the renormalon- 
based coefficients in a continuum scheme to the lattice 
scheme. Instead, one must map the perturbative coeffi¬ 
cients in the lattice scheme to a continuum scheme and 
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search for a scheme in which the mapped coefficients 
display a renormalon pattern. Once such a scheme is 
found one can then use Borel summation in that scheme 
to subtract the perturbative contribution to extract the 
nonperturbative power correction. 

In this paper we shall assume that such a scheme ex¬ 
ists and perform Borel summation using the bilocal ex¬ 
pansion of Borel transform |3l. To Borel-sum the di¬ 
vergent perturbative series accurately, one must have 
a precise description of the Borel transform in the do¬ 
main that contains the origin as well as the first renor¬ 
malon singularity in the Borel plane. The bilocal expan¬ 
sion, utilizing the known perturbative coefficients and 
the properties of the first renormalon singularity, recon¬ 
structs the Borel transform in the above domain of inter¬ 
est by interpolating the expansions about the origin and 
about the renormalon singularity. The Borel transform 
in the form of bilocal expansion implements the cor¬ 
rect nature of the first renormalon singularity as well as 
the perturbative coefficients that match the known co¬ 
efficients. The sum of the Borel-summed perturbative 
contribution and a dim-4 power correction is then fitted 
to the plaquette data. A good fit would then suggest the 
power correction be of dim-4 type. 

The Borel summation using bilocal expansion in the 
continuum scheme using the first A-loop perturbations 
of the plaquette is given in the form: 






/~*oo 

Ref 

F >"+ , 

Jo 



db, (12) 


where the integration is performed along the real axis 
on the upper-half plane. The essential idea of the bilo¬ 
cal expansion is to interpolate the two perturbative ex¬ 
pansions about the origin and about the renormalon sin¬ 
gularity to rebuild the Borel transform. Incorporating 
the renormalon singularity explicitly extends the appli¬ 
cable domain of the Borel transform from near the ori¬ 
gin to areas beyond the renormalon singularity. This 
scheme was used in summing the perturbative series 
for the static inter-quark potential as well as the heavy 
quark pole mass |lj,l5t]; The Borel-summed static poten¬ 
tial agrees very well with lattice results and the conver¬ 
gence in the pole mass case is extremely fast. 

N in (fTSI) denotes the normalization constant of the 
large order behavior and the coefficients hn are deter¬ 
mined so that the Borel transform reproduce the pertur¬ 
bative coefficients in the continuum scheme when ex¬ 
panded about the origin. Thus h„ depends on the con¬ 
tinuum perturbative coefficients as well as N. By def¬ 
inition, when expanded in 1//?, reproduces the 

perturbative coefhcients of the average plaquette to the 


A-loop order that were employed in building the Borel 
transform. The power correction can then be defined by: 

= (13) 

which, by definition, is free of perturbative coefficients 
to order N. 

Using the perturbation to 10-loop order of the plaque¬ 
tte we compute Pg°^(;S) in the continuum scheme de¬ 
fined by (lUi. The normalization N is treated as a fitting 
parameter, and in our scheme the parameters to be fitted 
are ri, r 2 , and N. 

Using the plaquette data for 6.0 < P < 6.8 from [@1 
and the relation between the lattice spacing a and p from 
static quark force simulation |f^ 

log(a/ro) = -1.6804- 1.733106-6) H- 

0.784906 - - 0.442806 - 6)^ (14) 

the fit gives Af = 165 and 

ri = 1.611, r2 = 0.246, (15) 

which values are substantially different from those in 
dTji. The result of the fit is shown in Fig. [1] which shows 
that the power correction is consistent with a dim-4 con¬ 
densate. The agreement improves as p increases, albeit 
with larger uncertainties; The deviation at low p (p < 6) 
may be attributed to a dim-6 condensate, which may be 
seen, though not presented here, by that adding a dim- 
6 power correction in the fit improves the agreement in 
the whole range of the plot. The error bars are from the 
uncertainty in the simulated perturbative coefficients of 
the plaquette. The uncertainty in the normalization con¬ 
stant does not appear to be large: for example, a varia¬ 
tion of 20% in N causes less than a quarter of those by 
the perturbative coefficients. 

From the fit we obtain a dim-4 power correction of 
Pnp ~ 1.6 {afrop. Because of the asymptotic nature of 
the perturbative series the power correction of the pla¬ 
quette is dependent on the subtraction scheme of the 
perturbative contribution, and thus our result may not 
be directly compared to those from other subtraction 
schemes. Nevertheless, it is still interesting to observe 
that the result is roughly consistent with 0.4 {alrof of 
iill and 0.7 {alrof of Our result turns out to be 
a little larger than those estimates; This may be partly 
accounted for by the fact that the existing results were 
from the fit in the low p range of y6 < 6, in which range 
the data are below our fitted curve. 

4. Comparison with latest result from 35-loop calcu¬ 
lation 

Recently Bali et al. Em computed the plaquette to 
35-loop order using stochastic perturbation theory and 
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Figure 1: log P^p vs. p. The solid line is for 4 logCo/ro) + 0.5. The plot shows 
the power correction should be of dim-4 type. 

observed a renormalon behavior of the coefficients at 
high orders, determining the normalization of the large 
order behavior. With the perturbative expansion in the 
asymptotic regime the power correction may as well be 
extracted by subtracting from the average plaquette the 
perturbative series U'uncated at the order the loop cor¬ 
rection becomes minimal. The power correction thus 
obtained was shown to be of dim-4 in accordance with 
the OPE ([TJ, and the gluon condensate was determined 
to be 

i—GG) ^ 0.077 GeV^ (16) 

n 

with an inU'insic uncertainty of 0.087 GeV"^. The nor¬ 
malization of the large order behavior in the lattice 
scheme was obtained as 

Mat = 4.2(+1.7)x 10^ (17) 

It is interesting to compare these numbers with our 
results. From the fit in the previous section we obtain 
the gluon condensate; 

i—GG) = ^ 0.14 GeV'^, (18) 

and the normalization N - 165+50 (30% uncertainty) 
in the continuum scheme corresponds to 

Mat = ye=^"'''M = 7.6(±2.3) x 10^, (19) 

which is, remarkably, consistent with (fTTl) . Comparing 
the gluon condensates (fTbl l and (fTSl l we see that they 
are in agreement within the intrinsic uncertainty. The 
gluon condensate is dependent on the prescription for 
the perturbative conU'ibution in the plaquette, and the 
difference between the two condensate values may be 
because they came from different prescriptions: Borel 
summation and truncated power series, respectively. 


5. Conclusions 

The renormalon subtraction procedure of 12] that led to 
a dim-2 condensate in the plaquette was reexamined. 
It is found that the continuum scheme employed in the 
procedure is far from a renormalon-dominated scheme 
and the procedure also fails a consistency check. As a 
consequence the power correction extracted is severely 
contaminated by perturbative contribution, to discredit 
the claimed dim-2 condensate. 

We then introduced a renormalon subtraction scheme 
that avoids the problems, in which the perturbative con¬ 
tribution is obtained by Borel-summing the perturba¬ 
tive series in a continuum scheme, employing the rebuilt 
Borel transform in the framework of the bilocal expan¬ 
sion. The power correction obtained in this procedure 
is of dim-4, in accordance with the OPE of the plaque¬ 
tte. The normalization of the large order behavior of the 
plaquette as well as the gluon condensate extracted are 
shown to be in agreement with the latest results from 
35-loop order calculations. 
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